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Abstract 



(N 

■ We study the issue of stability of static soliton-like solutions in 

<^ . some nondinear field theories which allow for knotted field configu- 

| rations. Concretely, we investigate the AFZ model, based on a La- 

grangian quartic in first derivatives with infinitely many conserved 
currents, for which infinitely many soliton solutions are known ana- 
<-j ', lytically. For this model we find that sectors with different (integer) 

topological charge (Hopf index) are not separated by an infinite energy 
barrier. Further, if variations which change the topological charge are 
allowed, then the static solutions are not even critical points of the 
energy functional. We also explain why soliton solutions can exist at 
all, in spite of these facts. In addition, we briefly discuss the Nicole 
model, which is based on a sigma-model type Lagrangian. For the 
Nicole model we find that different topological sectors are separated 
by an infinite energy barrier. 
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1 Introduction 



Nonlinear field theories, which allow for static soliton solutions of their field 
equations, have been studied intensively for many years due to their impor- 
tance in a variety of fields in theoretical and mathematical physics. Indeed, 
applications range from elementary particle theory, where they may serve 
as effective field theories providing particle-like solutions, to condensed mat- 
ter and solid state physics. One important criterion in the classification of 
these nonlinear field theories is given by the dimension of the base space 
(i.e., space-time) on which these fields exist. In the case of one space and 
one time dimension, static solutions have to solve a nonlinear ODE, which 
is usually much easier to solve than the nonlinear PDEs which result for 
solitons in higher dimensions. In addition, in 1 + 1 dimensions many rigor- 
ous results on nonlinear field theories and their soliton solutions are known, 
and a vast mathematical apparatus for the analysis of these models has been 
developed. Among those are inverse scattering, integrability, and the zero 
curvature representation, where the latter generalizes to field theory the Lax 
pair representation of integrable systems in + 1 dimensions (that is, systems 
with finitely many degrees of freedom). 

In higher dimensions (i.e., in d + 1 dimensional space-time for d > 1), on 
the other hand, much less is known about nonlinear field theories and their 
solutions. A generalization of the zero curvature representation to higher di- 
mensions was proposed in [1], and it was shown there that this proposal leads 
to nonlinear field theories with infinitely many conservation laws, realizing 
thereby the concept of integrability in higher dimensions in a concrete and 
well-defined manner. Besides this, some specific nonlinear field theories in 
higher dimensions have been studied with the intention of more phenomeno- 
logical applications. As far as static finite-energy solutions (solitons) are rela- 
vant in these applications, the selected nonlinear theories should, of course, 
really support such static solutions. One necessary condition for the existence 
of static finite-energy solutions is provided by the Derrick scaling argument 
(see, e.g., [2, 3]). The Derrick scaling argument simply says that static solu- 
tions with finite, nonzero energy cannot exist if it is possible for arbitrary field 
configurations with finite, nonzero energy to rescale the energy functional to 
arbitrarily small values by a base space scale transformation r — > fir, where 
r G IR d are the base space coordinates and //, with < /j, < oo, is the scale 
parameter. Let us be more precise for the class of theories we are interested 
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in. We assume that we have only scalar fields <£> a , a = 1 . . . N, which trans- 
form trivially under the above scale transformation, $ a — > $ a , and that the 
energy density for static configurations is a sum of terms each of which is 
homogeneous (of degree h) in space derivatives Vfc$ a . Then each such term 
£ h ($ , V fc $ ) contributes a term E h [$ a ] = J d d r£ h (<Z> a , V fc $ a ) to the energy 
functional, where EyJ^&a] transforms as ^[$ a ] — > fi d ~ h Ef l [^ a ] under the scale 
transformation. There are now several possibilities to obey the Derrick crite- 
rion for the possible existence of solutions. If the energy density only consists 
of one term, then necessarily d = h, that is, the degree of homogeneity in 
first derivatives, h, must be equal to the dimension of space, d. If the energy 
density consists of several terms, then at least one term must obey d — h>0, 
and at least one further term must have the opposite behaviour, d — h < 
(unless all terms obey the scale invariance condition d = h). All our concrete 
investigations will be for 3 + 1 dimensional space-time, therefore we assume 
d — 3 from now on. 

One well-known example of a field theory in d — 3 dimensions is the 
Skyrme model [4], which has applications in nuclear physics. It has the 
group SU(2) as field configuration space (target space), 1 and the fields are 
supposed to describe pions, whereas the soliton solutions are related to the 
nucleons. The energy functional of the Skyrme model consists of a sigma- 
model type term quadratic in first derivatives, with d — h = 3 — 2 = 1, and 
of a further term quartic in derivatives (the "Skyrme term") with d — h = 
3 — 4 = —1, so it obeys the Derrick criterion for the possibility of static 
solutions. Occasionally a third "potential" term is added to the Lagrangian 
density (and, consequently, to the energy density), which only depends on 
the fields, but not on derivatives. This term is usually assumed to make the 
pions massive. Soliton solutions have been found numerically for the Skyrme 
model both with and without the pion mass term [6, 7, 8]. 

Another well-known model with applications both to field theory and con- 
densed matter physics is the Faddeev-Niemi model ([9], [10]). This model 
may be derived from the Skyrme model by simply restricting its target space 
from SU(2) ~ S 3 to the two-sphere S 2 (e.g., the equator of the Skyrme 
S 3 , or SU(2)/U(1)). The target space of the Faddeev-Niemi model is two- 
dimensional, therefore the solitons are now line-like instead of point-like 
(where the soliton position may be defined, e.g., by the loci where the fields 

1 In particle physics, SU(3) has to be used, see e.g. [5]. 
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take values which are antipodal to their vacuum value). More precisely, the 
solitons of the Fadeev-Niemi model are links or knots. The corresponding 
topological index classifying a finite energy field configuration is a linking 
number (the Hopf index), whereas it is a winding number for the Skyrme 
model. The existence of soliton solutions in the Faddeev-Niemi model has 
been proven in [11], and confirmed by numerical calculations, e.g., in ([12] - 
[15]). 

The Faddeev-Niemi model has the Lagrangian density 



m 2 



C=—C 2 -XC 4 + C (1) 



where m is a constant with dimension of mass, A is a dimensionless coupling 
constant, £2 is 

(1 +uu) 

and £4 is 

(d^u d^u) 2 — (d^u dfju) (d v u d v u) 



Unit u r u 

^ = .,^ 2 , (2) 



'1+uuY 



(3) 



Further, u is a complex field which parametrizes the stereographic projection 
of the target S 2 , see Section 2 below for more details. Again we have allowed 
for a pure potential term £ (u,u). This term is not required by the Derrick 
criterion, but it may be necessary for some applications of the model (e.g., 
by providing a symmetry breaking of the target space symmetry group from 
SU(2) to U(l)). 

It is also possible to select Lagrangian densities which consist of one term 
only and are constructed from £4 or £ 2 exclusively, and obey the Derrick 
criterion h = 3. They are, however, necessarily non-polynomial. For £ 4 the 
appropriate choice is 

£afz = -(A) 1 - (4) 

This model has been introduced and studied by Aratyn, Ferreira and Zimer- 
man (AFZ) in [16], [17]. AFZ found infinitely many analytic soliton solutions 
for this model by using an ansatz with toroidal coordinates. The analysis of 
the AFZ model was carried further in ([18]), where, among other results, all 
the space-time and (geometric) target space symmetries of the AFZ model 
were determined, and, further, the use of the ansatz with toroidal coordinates 
was related to the conformal symmetry of the model (more precisely, of the 
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static equations of motion). It turns out that the AFZ model has infinitely 
many target space symmetries and, thus, infinitely many conservation laws. 
Moreover, it also shows classical integrability in a different sense, because 
the static field equations resulting from the ansatz with toroidal coordinates 
may be solved by simple integration. 

The other model, based solely on C 2 , is the Nicole model 

C m = (C 2 )l (5) 

This model has first been proposed by Nicole ([19]), and it was shown in 
the same paper that the simplest Hopf map with Hopf index 1 is a soliton 
solution for this model. (This model is, in fact, the restriction to S 2 target 
space of a non-polynomial £77(2) model which was studied first in [20] as a 
possible candidate for a pion model.) The Nicole model shares the confor- 
mal symmetry with the AFZ model and, therefore, the ansatz with toroidal 
coordinates may be used again to simplify the static field equations (to re- 
duce them to an ordinary differential equation). However, the Nicole model 
only has the obvious symmetries - the conformal base space symmetries (in 
the static case) and the modular target space symmetries, see [21]. Conse- 
quently, the field equations are no longer integrable, and the solutions are no 
longer available in closed, analytic form (except for the simplest case with 
Hopf index one 2 ). For a detailed investigation of soliton solutions within the 
ansatz in toroidal coordinates we refer to [23]. 

It is the main purpose of this paper to study stability issues of these two 
latter models with explicit solutions. Firstly, let us emphasize that the Der- 
rick criterion is a necessary condition for the existence of solutions, but by 
no means a sufficient one. So both the existence and the stability of soliton 
solutions have to be investigated independently. We find the more surprising 
results for the AFZ model. 3 In this model it turns out that there exists a sym- 
metry transformation (which maps solutions into solutions) which connects 
an arbitrary solution to the trivial vacuum solution u — 0. Consequently, 
soliton solutions with a nonzero topological index are not separated from the 
vacuum by an energy barrier. Further, soliton solutions are not even critical 
points of the energy functional, because there exists a nonzero first variation 

2 Exact solitons with higher Hopf index have been found in modified Nicole models [22] 
3 Warnings about its stability, based on the base space scale invariance r — > fir of the 
model, have been given in [3] . 
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of the energy functional into the direction of the symmetry transformation. 
This immediately poses an apparent paradox: (static) soliton solutions solve 
the variational equations of the energy functional, therefore they should be 
critical points of the energy functional by definition. If, on the other hand, 
the above-mentioned symmetry implies that there exists a non-flat direction 
of the energy functional for arbitrary fields, this seems to imply that finite en- 
ergy solutions cannot exist, much like the original Derrick scaling argument. 
Static soliton solutions of the AFZ model do exist, however, and have been 
constructed in ([16], [17]). In this paper, therefore, we not only show the ex- 
istence of the symmetry transformation connecting solutions to the vacuum, 
we also explain how the above-mentioned apparent paradox is resolved, by 
considering inherent source terms. 

In Section 2 we study the reduced energy functional which follows from an 
ansatz using toroidal coordinates. In Subsection 2.1 we perform this study 
for a simple toy model. The toy model has the same stability problems 
as the AFZ model, but much simpler field equations, therefore we mainly 
introduce it for illustrative purposes. This toy model may, however, be of 
some independent interest. It consists of a non-quadratic kinetic energy 
term, and such models have been studied in the context of "fake inflation" 
[24] (where the change in light propagation caused by the non-conventional 
kinetic term mimics a nontrivial geometry, as far as the propagation of light is 
concerned), and in the context of the so-called modified Newtonian dynamics 
(MOND), see, e.g., [25]. In addition, the toy model is equivalent to the 
electrostatic sector of a nonlinear theory of electrodynamics. This issue we 
discuss briefly in Subsection 3.1. In Subsection 2.2, at first we discuss the 
boundary conditions implied by the Hopf index for the ansatz within toroidal 
coordinates. Then we study the reduced energy functional of the AFZ model, 
using the analogy with the toy model. We establish the instability and show 
that the existence of solutions which are not critical points can be attributed 
to boundary term contributions to the reduced energy functional. In Section 
2.3 we briefly study the reduced energy functional of the Nicole model. We 
find that in this case different topological sectors are indeed separated by 
an infinite energy barrier, as one generally expects for topological soliton 
models. 

In Section 3 we study the full energy functional. In Subsection 3.1 we 
again investigate the toy model. We show that the apparent paradox is 
resolved by the presence of singular, delta-function like flux terms in a con- 
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servation equation. Further, we briefly introduce the above-mentioned non- 
linear electrodynamics. In Subsection 3.2 we investigate the AFZ model. It 
turns out that the apparent paradox related with the existence of solutions 
is resolved in this case by the existence of an inhomogeneous, singular flux 
term in just one of the infinitely many conservation equations of the model. 

Section 4 is devoted to the discussion of our results, emphasizing both the 
stability problems and the resolution of the apparent paradox. Further, we 
discuss possible generalizations of our results, and their application to and 
implications for further models. In Appendix A we calculate the Hopf index 
for a field within the ansatz in toroidal coordinates. 

2 The reduced energy functional 
2.1 A toy model 

Firstly, let us introduce a simple toy model which already shows the essential 
features we want to study in the sequel. We choose the energy functional 

E[$\ = J d 3 r(V$- V$)2 (6) 

where $ is a real scalar field. Here the non-integer power of | is chosen 
precisely as to make the energy functional invariant under a base space scale 
transformation r — > fir, avoiding thereby the usual Derrick scaling instabil- 
ity. There is, however, another scale transformation which does not leave 
invariant the energy functional. In fact, as both the energy functional and 
the resulting Euler-Lagrange equations are homogeneous in the scalar field 
<£>, the target space scale transformation $ — > A$ is a symmetry of the Euler- 
Lagrange equations which does not leave invariant the energy. As long as 
there exists no normalization condition for <3> which fixes the value of the scale 
factor A (and enters the energy functional, e.g., via a Lagrange multiplier), 
the same Derrick type scaling argument applies and seems to prevent the 
existence of finite energy solutions to the Euler-Lagrange equations. Never- 
theless, finite energy solutions exist, as we want to show now. We introduce 
toroidal coordinates (77, £, </?), rj e [0, 00], £, ip e [0, 2tt], via 

x = g _1 sinhr/cosy? , y = q~ x sinhr7sin(/3 

z = <7 _1 sin£ ; q = cosh 77 — cos £ (7) 
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in order to reduce the Euler-Lagrange equations to an ordinary differential 
equation (ODE). Further, we need the volume form 



dV = d 3 r = q 3 sinh rj dr] d£ dip (8) 

and the gradient 

V = (V77K + (V£)<% + (V</?H = + e& + — \^e v d v ) (9) 

where (e v , e$, e^) form an orthonormal frame in IR 3 . Assuming now that the 
scalar field <E> only depends on 77, the energy functional simplifies to 

roo 

£[$]=47r 2 / ^ sinh r/l^g 3 (10) 

J 

where & v = d n Q>. 

Remark: Due to the conformal base space symmetry of the Euler-Lagrange 
equations, these equations are compatible with the ansatz $ = $(77). It then 
follows from the principle of symmetric criticality that the resulting ODE is 
identical to the Euler-Lagrange equation which is derived from the reduced 
energy functional (10), see, e.g., [2, 3]. 

We prefer to introduce the new variable 

t = sinh 77 (11) 
which leads to the energy functional 

roo 

M$] = Aix 2 / dttil + t 2 )|$i| 3 . (12) 
Jo 

We now skip the absolute value signs by assuming that $t > 0V£, so that 
an infinitesimal variation of the scalar field leads to 

5E = E[<S> + 5$] - E[$\ 

roo 

~ 12tt 2 / dtt(l + t 2 )<& 2 t 5<$> t 
J 

= 12tt 2 



t(l+t 2 )^\ o - (13) 

roo rj 

12tt 2 / dt5<f>?-[t(l + t 2 )$ 2 } (14) 
Jo dt 



7 



and, therefore, to the Euler-Lagrange equation 



j t [t(i+m] = o as) 

with the solution 

$ t = C[t(l + t 2 )]-^ (16) 

where C is an integration constant. For C > it holds indeed that $ t > 0. 
The scalar field itself, 



m = cf , = = +c', (l?) 



can be expressed in terms of elliptic functions, but we do not need the explicit 
expression here. The energy of this field configuration is finite 

E = 4vr 2 C 3 [°° . dt = 32C 3 7rlr 2 (-) (18) 

(where T(-) is the Gamma function) and can take any positive value due to 
the arbitrary integration constant C. Further, for this field configuration the 
variation in the energy steming from the boundary term (13) is 

5E= 127r 2 C 2 [5$(oo) -5$(0)]. (19) 

Specifically, for a variation proportional to the field, 5$ = e$, this is nonzero, 

5E = 127r 2 C ,2 e[$(oo) - $(0)] = 96C ,2 e7rlr 2 (^). (20) 

In short, there exist finite energy solutions to the Euler-Lagrange equations in 
spite of the scaling instability under <3> — > A<3>, because such a variation 5$ = 
e$ produces a nonzero contribution to SE from the boundary term. It follows 
that the solution (16), despite being of finite energy, is not a critical point 
of the energy functional. The instability is in agreement with more general 
mathematical criteria for models without topological charges obtained by 
Rubakov and others, see, e.g., page 58 of Ref. [2]. 
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2.2 The AFZ model 



The degrees of freedom for both the AFZ and the Nicole model are given by 
a three- component unit vector field, that is, a map 

n(r) : H 3 -> S 2 , n 2 = 1 (21) 

where the tip of the unit vector field spans the unit two-sphere, or via stere- 
ographic projection 

1 ( , - ■( -\ i -\ ni +m 2 , . 

n = : — [u + u, —i{u — u), 1 — -u-uj ; w = . (22) 

1+ I u \ 2 1 + n 3 

by a complex field 

u(r) : IR 3 -> <Cq (23) 

where our conventions are such that the projection is from the south pole 
to the equatorial plane of the two-sphere. Further, (D is the one-point com- 
pactified complex plane. 

Topological solitons, i.e., static finite energy solutions with a non-zero 
integer value of the corresponding topological charge (the Hopf index) have 
been found for both models, so let us briefly describe the corresponding 
topological map (Hopf map). A Hopf map is a map : S* 3 — > S 2 or, via 
the stereographic projection from S 3 to one-point compactified IRq, a map 
(f) : IRq — > S 2 , and is characterized by the integer Hopf index (for details we 
refer to the Appendix A of Ref. [23]). Therefore, the above fields n and u 
can be identified with Hopf maps provided that the following two conditions 
hold. Firstly, the fields have to obey 

lim n = n = const , lim u = u = const (24) 

r|^oo |r|^oo 

so that they are, in fact, defined on one-point compactified Euclidean three- 
dimensional space IRq. Secondly, the fields have to take values in the full 
target spaces S 2 and (D , respectively, in such a way that the preimage of 
each point in target space is a closed line in the base space IRq. 

Both the AFZ and the Nicole model have a conformally invariant energy 
functional and conformally symmetric Euler-Lagrange equations, therefore 
the Euler-Lagrange equations are, again, compatible with the ansatz 

u = f(ny* +imp (25) 
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using toroidal coordinates. The limit |r| — > oo corresponds to the limits £ — > 
and ?7 — > 0, therefore / has to obey lim^o/ (v) = or lim^o f(v) = 00 
in order that u be defined on IRq. As u and 1/u are related by a symmetry 
transformation in both models, we may asssume for a true Hopf map that 
f(r] — 0) = without loss of generality. For u to take values in the whole 
target space (Do now requires that / takes values on the whole positive real 
semiaxis IRq . Further, / has to take the value oo at r\ = oo, because the 
pre-image of r\ = oo is a circle (the unit circle in the (x, y)-plane), whereas 
the preimages of finite i] = const, are surfaces (tori). So if we want the 
preimage of the point u = oo to be a circle, u has to take the value oo at 
rj = oo. In short, for a u within the ansatz (25) to be a genuine Hopf map, 
/ has to obey 

/(77 = 0) = 0, f(r ] = oc) = oc or /(* = 0) = , f(t = oo) = oo 

(26) 

where t = sinh 77, as above. In Appendix A we derive the same result from 
an analytic expression for the Hopf index by proving that integer valuedness 
of the Hopf index precisely requires the boundary conditions (26). 
The energy functional for the AFZ model is 

and, for the ansatz (25) and using symmetric criticality, 

£afz[/]=4tt dr)smhr)\ (I + py — J ' ^ > 

For reasons of simplicity we now restrict to the case m = n (the discussion 
for n 7^ m is completely analogous). Further we define F = f 2 and find for 
the energy functional 

3 

The solution of the Euler-Lagrange equation which obeys the boundary con- 
dition F(t = 0) = 0, F(t = 00) = 00 required for a Hopf map is F = t 2 , see 
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below. However, field configurations which deviate from these boundary con- 
ditions give perfectly valid, nonsingular energy densities with finite energy 
E.g., for F = c + t 2 , c > 0, the energy is 

E AFZ (c ) = An 2 ml f°° dtt(l + t 2 )(l + c + t 2 )' 3 = 7r 2 (2m)§ 2 ± \ . (30) 
Jo [1 + CqY 

Exactly the same result is obtained for F = t 2 /(l + c t 2 ), which deviates 
from the boundary condition for a Hopf map at t — oo. 

Next, we want to find the Euler-Lagrange equation. We introduce the 
function 



which leads to the energy functional 



3 fOO I + f 3 

E A fz[G] =Aix 2 m* / dt-^\G t \* (32) 

which is quite similar to the energy functional of the toy model in the pre- 
ceding section. As the energy functional is homogeneous in G, we again have 
that the target space scale transformation 

G -> AG (33) 

is a symmetry of the Euler-Lagrange equation which does not leave invariant 
the energy. 

Remark: as long as we only consider the energy functional (32) by itself, 
A may take arbitrary values. However, the relation (31) requires, for positive 
semidefinite F, that — G G [0, 1] and, therefore, restricts the possible values 
of A to < A < 1. 

For the variation of the energy functional we easily find (assuming again 
G t >0Vt) 



SE AFZ = E AFZ [G + 5G\ — E AFZ [G] 
~l + t 2 



2 2 
= 07T ffl2 



Vt 



G?5G 



(34) 



67r 2 m§ / dt8G%-Y-^Gl\ (35) 
Jo dt [ Ji J V ; 
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and therefore the Euler-Lagrange equation 



(39) 



d 1 + t 2 1 

s^n- < 36 > 

with the solution 

G ' = 1JTTW (37) 

G = -JL^-D (38) 

where C and D are non-negative real integration constants. Further, 

l-C-D + t 2 (l-D) 
C + D(l + t 2 ) 

is positive semidefinite by construction, therefore C and -D are restricted to 

C + D<1 A D<1. (40) 

For C = 1, D = we find the solution F = t 2 which obeys the boundary 
conditions for a Hopf map. For other values of C and D we find solutions 
F which do not obey these boundary conditions and, therefore, do not give 
rise to Hopf maps. Nevertheless, they are regular configurations which have 
finite energy. Specifically, for C strictly less than one the energy is lowered, 

E AFZ = ^-K 2 mic'i. (41) 

Again, for variations proportional to the solution, 5G = eG, the boundary 
term (34) gives a nonzero contribution to the variation of the energy, 

SE AFZ = 3V2ir 2 mhC%. (42) 

Therefore, the solution to the Euler-Lagrange equation is again not a critical 
point of the energy functional. 

At this point we may ask to which symmetry transformation corresponds 
the target space scale transformation G — * AG in terms of the original fields u 
and u. This transformation is, in fact, already well-known (the infinitesimal 
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version has been found in [18], and the finite version has been calculated in 
[26]), and it does not depend on the specific choice of toroidal coordinates, but 
rather is a pure target space transformation. It is based on the observation 
that the action density £ 4 is just the square of the pullback (under the map 
u) of the area twoform on the target space S 2 , 

„ . dudu , 
(1 + uu) 2 

The transformation we search is the transformation u — > v(u,u) such that 
the area twoform is mapped to 

dudu dvdv . 9 dudu , , , . 

= A (44) 



(1 + uu) 2 (1+vv) 2 (1+uu) 2 ' 

If we introduce the real coordinates on target space u = F x l 2 e ia (angle and ra- 
dius squared on the Euclidean plane; here we do not assume a specific variable 
dependence of the modulus squared F) and assume that v = (F) 1 / 2 (F)e U7 
(i.e. u and v have the same argument, and the transformation only affects 
the modulus) then we get the equation 

F'(F)dFda t2 dFda 

= A (45) 



or 



with the solution 



(1 + F) 2 (1 + F) 2 

F r A 2 
(1 + F) 2 ~ (1 + F) 2 

1 A 2 



(46) 



r, + C ( 47 ) 

l + 1 + F v ; 

where c is a constant of integration. For c = we just get G = A 2 G, that 
is, the scale transformation (33) with A = A 2 . If, instead, we impose the 
boundary condition that v = for u = then c = 1 — A 2 and the solution is 

A 2 F , . 

(48) 



A 2 + (1 + F)(l — A 2 ) 

or 



Am , 

(49) 



[A 2 + (1 + mm)(1-A 2 )]V 
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Observe that A, again, is restricted to A < 1 if u is a genuine Hopf map 
which covers the whole target space. Further, v is no longer a true Hopf map 
for A < 1. 

Remark: the transformation (49) is a target space symmetry transfor- 
mation of the field equations, i.e., a target space transformation which maps 
solutions u into solutions v. It is, however, not a symmetry of the action, but, 
instead, rescales the action density £ 4 by £ 4 — > A 4 £ 4 . It is, therefore, not 
a Noether symmetry and does not define a Noether current and conserved 
charge. On the other hand, a combination of the target space transformation 
(49) and the base space scale transformation (dilatation) (t, r) — > (fit, fir), 
fi = A -4 , leaves the action invariant and, therefore, gives rise to a Noether 
current and conserved charge, see [18] for details. 

The transformation (49) is a pure target space transformation and, there- 
fore, the rescaling of the action density £ 4 by £ 4 — > A 4 £ 4 is independent of 
the base space. As a consequence, stability problems which are similar to 
the ones discussed here can be expected for all lagrangian densities which are 
homogeneous in £ 4 , for arbitrary base spaces. This issue will be discussed 
further in the last section. 

2.3 The Nicole model 

Here we just want to demonstrate briefly that, at least for field configurations 
belonging to the ansatz (25), soliton solutions of the Nicole model with non- 
zero Hopf index are indeed separated by an infinite energy barrier from the 
trivial sector, and that field configurations which do not obey the boundary 
conditions necessary for an integer Hopf index have infinite energy. The 
energy functional of the Nicole model is just a non-integer power of the CP 1 
model Lagrangian, 

For the ansatz (25), and for m — n, which again we assume for simplicity, 
this reduces to 

3 

roo ( f 2 \ 2 poo 

Em[f]=J Q dtt(l + t 2 )(l + f)- 3 [f? + rn 2 ^) = J q dt£(t). (51) 
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Now let us assume that /(0) = c > 0. For the density S(t) this implies 

3 

]im£(t) ~ limt(l + c^)- 3 (/ t 2 + m 2 |l ' > ^~ 2 7T~~~|\s (52) 

*-o w *-o v uy t 2 y - (l + c 2 ,) 3 v 7 

which has a nonintegrable singularity at £ = 0. Equivalently, in the limit 

t — > oo we have, assuming /(oo) = Coo, 

3 

Urn £(*) = «•(! + cD- (/? + m^) 5 > (53) 

Therefore, the density is bounded from below by a nonzero constant, and its 
integral again will be infinite. 

It follows that for field configurations which obey boundary conditions 
such that the corresponding Hopf index is non-integer, the energy in the 
Nicole model automatically is infinite. In other words, it is a necessary 
condition for a finite energy field configuration to have integer Hopf index. 
The condition is not sufficient, however. For fields with the leading behaviour 
\im t ^o£(t) — t a °, linit-nxj £ (t) ~ t a °° , for instance, it may be checked easily 
that finite energy requires > 1/3 and > 1/3, which is more restrictive 
than just /(0) = and /(oo) = oo. 

3 Full energy functional and nonzero fluxes 
3.1 Toy model 

The variation of the energy functional of the toy model for an arbitrary scalar 
field is 

5E = E[Q + 6$] - E[<S>] 

= 3 J d 3 rV • (5$|V$|V$) (54) 

- 3 j d 3 r5$V • (|V$|V$) (55) 

and the resulting Euler-Lagrange equation is just the conservation equation 

V • (|V$|V$) = V • J = 0. (56) 
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The nonzero contribution of the solution (17) to the boundary terms of 
the variation SE of the reduced energy functional indicates the presence of 
nonzero delta-function like source terms in the conservation equation (56). 
Further, these source terms we expect to be concentrated at the loci of the 
boundaries of the reduced system, that is, at rj = (the z axis) and at rj = oo 
( the unit circle in the (x, y) -plane) . This we want to investigate in the sequel. 
For the solution $(77) of the last section with = C sinh^ 5 ^ the current J 
is 

J -^-C3^ (5?) 

We now want to calculate the flux of this current through a torus rj = const. 
The surface element of a torus rj = const, is just 

Slllll Tj 

dT, T = — - — e v d^dip (58) 
q l 

therefore the flux simply is 

//•27T /•27T 
d£ T • J = / dU dipC 2 = 4ir 2 C 2 . (59) 
Jo Jo 

So there is some nonzero total flux emerging from the z axis and streaming 
towards the unit circle 

C = {x G IR 3 : z = A r 2 = 1}. (60) 

The total flux escaping to infinity is zero, as may be checked by explicit 
calculation. However, this already follows from the fact that the total charges 
distributed along the z axis and along the unit circle C are equal in magnitude 
and opposite in sign. The line charge density along the circle C is constant, 
as is obvious for symmetry reasons. The line charge density along the z axis 
may be calculated by calculating the flux through a cylinder around the z 
axis with infinitesimally small radius. For infinitesimally small radius, the 
top and bottom of the cylinder do not contribute. The surface element of 
the cylinder mantle in cylinder coordinates (p, ip, z), p 2 = x 2 + y 2 , is 

dT,M = pipdfdz (61) 
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Further, in the limit p — > 0, e v ~ e p , and therefore 

/rZ2 r2ir 
dT,M ■ J = lim C 2 dz dtpp 
p^o J Z1 Jo 



sinh?7 



= 2ttC 2 



+ z l 



where we used 



2 sinh 77 

P = 



and 



sinh 2 rj 



4p 2 



4^ 2 + (p 2 + - 1) 



sinh 77 
hm 

p^O p 



1 + z 2 



(62) 



(63) 



(64) 



Consequently, the solution $(77) solves in fact the inhomogeneous equation 
V • J = 4nC 2 Q = 4nC 2 (s(x)S(y)^^ - 5(z)5(p 2 - 1)) , (65) 
and the total charge of the density distribution Q is zero, 

J dhAitC 2 Q = J d 3 rV • J = 0, (66) 

as may be checked easily. 

Remark: The observation that the toy model is equivalent to the purely 
electrostatic sector of a nonlinear theory of electrodynamics may add some 
additional interest to it. Indeed, take the Lagrangian density 



1 ( F nF a P \ 2 

r---F F*™ I p 

4 V 2ft 4 



(67) 



where ft is a constant with the dimension of mass, and F^ u = d li A v — d u A^, 
= (V, A). This leads to the field equation (Gauss law) 



F 



F n aF a ^ 



2ft 4 







(68) 



or, in the electrostatic case (A = 0, V = V(r), E = — VV), to 

V-(|W|W) = (69) 
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which is identical to Eq. (56) when the potential V is identified with the 
scalar field $ of the toy model. Therefore, we have the electrostatic solutions 

Vt = i ° (70) 

and 

-> C 

E = -qV v e v = -qVl + t 2 V t e v = -q^=e v . (71) 



The electric induction in this model 



D-^ 2 \E\E (72) 

is proportional to the current J of the toy model and, therefore, obeys the 
inhomogeneous divergence equation 

of 2 

V-5=-— Q (73) 
where the singular line charge density Q is defined in Eq. (65). 

3.2 AFZ model 

The variation of the energy functional of the AFZ model is 
SEafz = E[u, u + 5u] — E[u, u] 

= - J d 3 rV ■ (5u(l + uu)- 2 L) (74) 



J d 3 r5u(l + uu)~ 2 V L (75) 



where 



L = (l + uu^H-^K (76) 

H = (Vu • Vm) 2 - (Vm) 2 (Vm) 2 (77) 

K = (Vu • W)Vw - (Vm) 2 Vm. (78) 

The resulting Euler-Lagrange equation is the conservation equation 

V • L = 0. (79) 
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We want to investigate whether, again, there are delta-function like source 
terms present at the r.h.s. of this conservation equation. However, for the 
AFZ model things are complicated by the fact that, due to the infinitely 
many symmetries of the model, there exist in fact infinitely many currents 
such that their conservation is equivalent to the Euler-Lagrange equations. 
Indeed, with the help of the identity 

V-u • K = (80) 

it easily follows that together with the conservation equation (79) we have 
the infinitely many conservation equations 

V • L c = V • CL = (81) 

where ( = ((u) is an arbitrary function of u only. The currents are, 

in general, complex, but it is easy to construct an equivalent set of real 
conserved currents, 

l g = -i (g u L - g u L) (82) 

where Q u = d u Q, etc., and Q = Q(u, u) is an arbitrary real function of u and 
u. The real currents L e are just the Noether currents of the area-preserving 
target space diffeomorphisms, which are wellknown symmetries of the AFZ 
model [18, 27]. 

For the ansatz (25) for u in toroidal coordinates (and m = n), the current 
(76) is 



3 1 

a 2 



tanhr] J ytanhr? J \ sinhr^ 

(83) 

This current will give a zero total flux through a torus r\ = const, because of 
the angular factor exp[im(£ + </?)]. By inspection, it is obvious that a nonzero 
flux will exist for the current with ( = u" 1 , because then the angular 
factor is absent. Multiplying by the surface element of a torus r\ = const, we 
find 



- \ 2 



t i jri 3 COSll2 ( ff v 

L« • Gt£ T = 2m 2 



sinh^ V(l + / 2 ) 2 
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= V2: 



m 2 • 



1 + F t 
Vt {(1 + F)" 



n- - 



where we used the solution 

F t 



(1 + F)' 



= G t 



C t 
2"(l + t 2 ) 2 ' 



(84) 



(85) 



and therefore the flux +n 2 mh\fC through an arbitrary torus of rj = const. 

Remark: the current with nonzero flux may be equally found among the 
real currents L e for the choice 



G = Q { = -(hi if — \nu). 



(86) 



The current L^ f is identical to the current J of the toy model and, therefore, 
obeys the inhomogeneous divergence equation 



V ■ L gl = AnC^ml (5(x)5(y)—^ - 5(z)5(p 2 - 1)) . 



(87) 



4 Discussion 



The main result of this paper is the observation that in the AFZ model for 
any solution u there exists a one-parameter family of solutions v(u, A) which 
connects the given solution u at A = 1 to the vacuum solution v(u, 0) = 0, see 
Eq. (49). Therefore, sectors of field configurations with different topological 
(Hopf) index are not separated by an infinite energy barrier, and there exist 
well-defined (i.e. single-valued, nonsingular) field configurations with finite 
energy in the model which do not correspond to integer Hopf index. For 
static solutions (i.e., solitons), this observation immediately implies some 
problems. The first problem is that the mere existence of soliton solutions 
under these conditions poses an apparent paradox. Indeed, the static field 
equations which the solitons obey are just the Euler-Lagrange equations for 
the variational problem of the energy functional Eafz[u, u] and, therefore, if 
u is a soliton solution, it should hold that 



d_ 

de 



E AF z[uo,u + e5u] 



6=0 



d_ 

de 



Eafz[uq + e5u,u } 



= 



(88) 



6=0 
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for arbitrary 5u. But in our case this is not true. We have that E A fz[v, v] 
A 3 Eafz[u, u] and therefore 



+ 

e=0 



-^-E AFZ [u + e5u,u } 
de 



= 3E AFZ [uo,u ] (89) 



6=0 



for A = 1 — e and 5u = u (l + u u ). The resolution of this apparent para- 
dox can be understood in different ways. For the reduced energy functional 
within the ansatz in toroidal coordinates (see Section 2), it turned out that 
in addition to the bulk contribution to the variation of the energy functional 
(which is zero for a solution) there are boundary contributions, which are 
nonzero for certain variations, even when the variations are performed about 
a soliton solution. 

For the full energy functionals the nonzero contributions to the variation 
about a solution cannot be attributed to boundary terms, because there are 
no boundaries in the base space IR 3 , and the soliton solutions are behaving 
well in the limit r — > oo. For the toy model, the resolution of the paradox is 
quite simple. The static field equation corresponds to a conservation equa- 
tion V • J = for some current J, whereas the "solution" obeys, in fact, an 
inhomogeneous conservation equation instead, with some delta-function like 
sources and sinks on the r.h.s. For the AFZ model, the resolution is slightly 
more subtle. The static field equation is again equivalent to a conservation 
equation V • L = 0, but now the well-known soliton solutions do not in- 
duce an inhomogeneous source term in this conservation equation. However, 
due to the infinitely many symmetries of the AFZ model, there formally 
exist infinitely many more conservation equations V • = for soliton solu- 
tions. It turns out that, for nontrivial soliton solutions, one of these infinitely 
many conservation equations is, in fact, inhomogeneous and contains, again, 
a singular flux term on the r.h.s., see Eq. (87). Further, the strength of 
the delta-function like inhomogeneous term is not varied in the variational 
derivation of the field equations. Therefore, the value of A, which plays the 
role of an integration constant for solutions within the ansatz in toroidal co- 
ordinates (A = y/C, see Eqs. (85) and (87)), remains fixed, which explains 
the existence of solutions to the variational equations. 

The next issue to be discussed is, of course, the stability of solutions (that 
is, the stability under small perturbations). Naivly, one might assume that 
the presence of the inhomogeneous term in one of the conservation equations 
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stabilizes the solutions by fixing the value of A. However, this might not be 
correct for the following reasons. Firstly, the inhomogeneous term does not 
correspond to a nonzero charge, because the spatial integral of the inhomo- 
geneity is zero, see Eq. (66). Secondly, the fact that the inhomogeneous term 
is not varied by the variation of the energy functional does not imply that 
it is conserved under dynamical evolution of the system. The conservation 
under dynamical evolution would require that A (or a related quantity, like, 
e.g., the Hopf charge, which varies under the transformation (49)) may be 
expressed by one or more of the conserved charges of the theory. For static 
solutions this most likely does not happen, however, because most of the 
conserved charges are trivially zero for static solutions (e.g., all the infinitely 
many conserved charges related to the infinitely many area-preserving target 
space diffeomorphisms contain time derivatives and are, therefore, zero for 
static solutions). We conclude that the stability of the soliton solutions in 
the AFZ model is problematic. 

Although the explicit calculations have been done for the static AFZ 
model with Lagrangian (£4)4 and base space IR 3 , some conclusions may be 
drawn for more general models based on the Lagrangian £ 4 and with different 
base spaces. For instance, static soliton solutions have been found in [28] for 
the Lagrangian £4 and the three-sphere S 3 as base space (that is, space-time 
S* 3 x IR). These static solutions solve the Euler-Lagrange equations which 
follow from varying the static energy functional 



(where dVs3 is the volume element on S* 3 ) and have finite energy, therefore the 
same arguments as above apply. Due to the symmetry (49) these solutions 
are connected to the vacuum and, in addition, are not critical points of the 
energy functional. Therefore, the mere existence of these solutions requires 
the presence of inhomogeneous singular flux terms in some (at least one) of 
the infinitely many conservation equations. This line of arguments remains 
valid for static finite energy solutions for Lagrangians which are arbitrary 
powers of the Lagrangian £4, and for arbitrary base spaces. 

The situation is slightly different for time-dependent solutions. Indeed, 
time dependent solutions for the Lagrangian £4 have been found both for base 
space-time IR 4 and S 3 x IR in [29] and [30], respectively. These solutions have 
integer Hopf charge and are particle-like, that is, they have finite energy and 




(90) 
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infinite action. These solutions are again connected to the vacuum by the 
symmetry transformation (49), but the apparent paradox discussed above 
is absent for these time dependent solutions. The energy transforms like 
E — > A A E under the transformation (49), but the time dependent field equa- 
tions are not derived by varying the energy functional, and time dependent 
solutions are not required to be critical points of the energy functional. The 
action / transforms like the energy under the symmetry transformation (49), 
I — > A 4 /, so the same apparent paradox would exist for time dependent finite 
action solutions ("instantons"). But it does not exist for time dependent so- 
lutions with infinite action (particle-like solutions). Therefore, the existence 
of inhomogeneous singular flux terms in some conservation equations is not 
required by the existence of these time dependent solutions (although we 
do not know whether such inhomogeneous singular flux terms do or do not 
show up for the time dependent solutions). Even the issue of stability may 
be less problematic for these time dependent solutions. The point is that the 
infinitely many conserved Noether charges, which exist in the model as a con- 
sequence of its infinitely many symmetries, are nontrivial (i.e., nonzero) for 
time dependent solutions. It is quite possible that these infinitely many con- 
servation laws stabilize the time dependent solutions in spite of the fact that 
those solutions are connected to the vacuum by a symmetry transformation. 

Another interesting feature of models based on the Lagrangian density 
£4, which is related to the existence of finite energy configurations with non- 
integer Hopf index, is the rich vacuum structure of such models. Indeed, a 
complex field u(h) which only depends on one real function h(t, r) gives zero 
when inserted into the Lagrangian density £4, implying that the vacuum 
manifold is infinite dimensional. This may be understood in geometric terms 
by the observation that the quartic Lagrangian £ 4 (and the related energy 
density) is just the square of the pullback (under the map u) of the area two- 
form on target space S* 2 . A field u(h) of the above type is a map IR 3 — > M 1 — > 
S 2 , where Ai 1 is a one-dimensional manifold, and the pullback of a two- form 
onto a one-dimensional manifold is zero. As a consequence, static fields u 
which have a functional dependence u(h) in the limit r — > 00 will have finite 
energy. In the AFZ model and within the ansatz of toroidal coordinates, for 
instance, it happens that lim r _ >0O 'u = c exp(irrup) whenever f(rj) obeys the 
boundary condition /(0) = c . 

The rigorous results, as well as the detailed arguments summarized in 
this section for some nonlinear models in 4 dimensions, are the main con- 
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tribution of our work. These models are special, but they are related to 
and share many properties of some relevant fundamental and effective field 
theories. The Lagrangian £ 4 , for instance, is just the quartic part of the 
Faddeev-Niemi model. The stability of this term has been investigated re- 
cently in [31], together with some generalizations (Lagrangians based on the 
Kaehler two-form on various target spaces). However, the stability analysis 
of the quartic Lagrangian £ 4 in [31] focused on the quadratic variation of 
the energy functional, whereas we found that already the linear variation is 
problematic. Our general conclusion, therefore, is that the stability of the 
extended solutions, which is often taken more or less for granted, is an ex- 
tremely subtle and difficult issue which, as shown here, can be understood 
in detail by a careful analysis in some cases. Our methods can certainly be 
useful for that purpose. Note also that in numerical analysis the restriction 
to a fixed topological sector is often imposed from the outset by assuming 
the corresponding boundary conditions. Under these conditions, it would not 
be possible via numerical analysis to detect stability problems of the type 
discussed here. 



Appendix A 

Here we want to demonstrate briefly from the analytic expression for the 
Hopf index that for fields u within the toroidal ansatz (25) the condition 
that u has integer Hopf index precisely requires that / obeys the boundary 
condition (26). The analytic expression for the Hopf index is 

Q =iLJ d3rAe (9i » 

where B is the Hopf curvature 

g = 2 V«xV« 

i(l + uu) 2 v ; 

and A is the gauge potential for the "magnetic field" B, B = V x A. There 
is no local expression for A in terms of the Hopf map u alone, but when u is 
expressed with the help of a four component unit vector field e a , a — 1, ... ,4 
like 

u = , (93) 

e 3 + ze 4 
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then an explicit, local expression for the gauge potential A in terms of the 
e a exists and is given by 



A = -[(ei - ie 2 )V(ei + ie 2 ) + (e 3 - ie 4 )V(e 3 + ie 4 )}. 
i 



In our case we may choose 



ei + ie 2 



./ 



imip 



e 3 + ie 4 = 



1 



and, therefore, 



A 2 1 

A = -q 



( vu 

B = 4q 



e v + ~ 



imf 2 



sinh?7 



e v — me^ 



2 ffy 
(1 + / 



m 



2\2 



and 



Q = J dr]d£d(p8mn— 



ffr. 



sinhr/ 



mn 



1 + F(rj) 



J o 



(94) 
(95) 

(96) 
(97) 

(98) 



+ P) 2 

(remember F = f 2 ). It follows that the conditions for integer, nonzero 
Hopf index are precisely /(0) = and /(oo) = oo (or the inverse conditions 
/(0) = oo and /(oo) = 0; however, as u — > (1/u) is a symmetry transforma- 
tion for both models, these cases are equivalent). 
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